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Abstract 



Aims. We investigate whether the morphology of bilateral supernova remnants (BSNRs) observed in the radio band is determined mainly either 
by a non-uniform interstellar medium (ISM) or by a non-uniform ambient magnetic field. 

Methods. We perform 3-D MHD simulations of a spherical SNR shock propagating through a magnetized ISM. Two cases of shock propagation 
are considered: 1) through a gradient of ambient density with a uniform ambient magnetic field; 2) through a homogeneous medium with a 
gradient of ambient magnetic field strength. From the simulations, we synthesize the synchrotron radio emission, making different assumptions 
about the details of acceleration and injection of relativistic electrons. 

Results. We find that asymmetric BSNRs are produced if the line-of-sight is not aligned with the gradient of ambient plasma density or with 
the gradient of ambient magnetic field strength. We derive useful parameters to quantify the degree of asymmetry of the remnants that may 
provide a powerful diagnostic of the microphysics of strong shock waves through the comparison between models and observations. 
Conclusions. BSNRs with two radio limbs of different brightness can be explained if a gradient of ambient density or, most likely, of ambient 
magnetic field strength is perpendicular to the radio limbs. BSNRs with converging similar radio arcs can be explained if the gradient runs 
between the two arcs. 

Key words. Magnetohydrodynamics (MHD) - Shock waves - ISM: supernova remnants - ISM: magnetic fields - Radio continuum: ISM 



1. Introduction 

It is widely accepted that the structure and the chemical abun- 
dances of the interstellar medium (ISM) are strongly influ- 
enced by supernova (SN) explosions and by their remnants 
(SNRs). However, the details of the interaction between SNR 
shock fronts and ISM depend, in principle, on many factors, 
among which the multiple-phase structure of the medium, its 
density and temperature, the intensity and direction of the am- 
bient magnetic fields. These factors are not easily determined 
and this somewhat hampers our detailed understanding of the 
complex ISM. 

The bilateral supernova re mnants (BSNRs, Gaenslerlll998 ; 



1987, 



also call ed "barrel-shaped," Kesteven & Caswell 
"bipolar". lFulbright & Reynolds! 199ol) are considered a bench- 
mark for the study of large scale SNR-ISM interactions, since 
no small scale effect like encounters with ISM clouds seems 
to be relevant. The BSNRs are characterized by two opposed 
radio-bright limbs separated by a region of low surface bright- 
ness. In general, the remnants appear asymmetric, distorted and 
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elongated with respect to the shape and surface brightness of 
the two opposed limbs. In most (but not all) of the BSNRs 
the symmetry axis is parallel to the galactic plane, and this has 
been interpreted as a difficulty for "intrinsic" models, e.g. mod- 
els based on SN jets, rather than for "extrinsic" models, e.g. 
models based o n properties of the surrounding galactic medium 
dGaenslerlll998l) . 



In spite of the interest around BSNRs, a satisfactory and 
complete model which explains the observed morphology and 
the origin of the asymmetries does not exist. The galactic 
medium is supposed to be stratified along the lines of con- 
stant galactic latitude, and characterized by a large-scale am- 
bient magnetic field with field lines probably mostly aligned 
with the galactic plane. The magnetic field plays a three-fold 
role: first, a magnetic tension and a gradient of the magnetic 
field strength is present where the field is perpendicular to the 
shock velocity leading to a compression of the plasma; sec- 
ond, cosmic ray acceleration is most r apid where the field lines 
are pe rpendicular to the shock speed dJokipiill 19871 lOstrowski 
1988b ; third, the electron injection could be fa vored where 



the magnet ic fie l d is p arallel to the shock speed (Elli son et al 



1995b . iGaenslerl (11998b notes that magnetic models (i.e. those 
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considering uniform ISM and ordered magnetic field) cannot 
explain the asymmetric morphology of most BSNRs, and in- 
vokes a dynamical model based on pre-existing ISM inhomo- 
geneities, e.g. large-scale density gradients, tunnels, cavities. 
Unfortunately, the predictions of these ad-hoc models have 
consisted so far of a qualitative estimate of the BSNRs mor- 
phology, with no real estimates of the ISM density interacting 
with the shock. Moreover, most likely also non-uniform ambi- 
ent magnetic fields may cause asymmetries in BSNRs, without 
the need to assume ad-hoc density ISM structures. Two main 
aspects of the nature of BSNRs, therefore, remain unexplored: 
how and under which physical conditions do the asymmetries 
originate in BSNRs? What is more effective in determining the 
morphology and the asymmetries of this class of SNRs, the 
ambient magnetic field or the non-uniform ISM? 

Answering such questions at an adequate level requires 
detailed physical modeling, high-level numerical implementa- 
tions and extensive simulations. Our purpose here is to inves- 
tigate whether the morphology of BSNR observed in the ra- 
dio band could be mainly determined by the propagation of 
the shock through a non-uniform ISM or, rather, across a non- 
uniform ambient magnetic field. To this end, we model the 
propagation of a shock generated by an SN explosion in the 
magnetized non-uniform ISM with detailed numerical MHD 
simulations, considering two complementary cases of shock 
propagation: 1) through a gradient of ambient density with a 
uniform ambient magnetic field; 2) through a homogeneous 
isothermal medium with a gradient of ambient magnetic field 
strength. 

In Sect. [2] we describe the MHD model, the numerical 
setup, and the synthesis of synchrotron emission; in Sect. [3] we 
analyze the effects the environment has on the radio emission 
of the remnant; finally in Sect. H] and [5] we discuss the results 
and draw our conclusions. 



2. Model 

2.1. Magnetohydrodynamic modeling 

We model the propagation of an SN shock front through a 
magnetized ambient medium. The model includes no radia- 
tive cooling, no thermal conduction, no eventual magnetic field 
amplification and no effects on shock dynamics due to back- 
reaction of accelerated cosmic rays. The shock propagation 
is modeled by solving numerically the time-dependent ideal 
MHD equations of mass, momentum, and energy conservation 
in a 3-D cartesian coordinate system (x,y,z): 



% + V • (pu) = , 
ot 

+ V • (puu - BB) + VP, = 

ot 



dpE 
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where 



B 2 

P«=P+ — , 
2 



1 , 1 |B| 2 

2 2 p 



are the total pressure (thermal pressure, P, and magnetic pres- 
sure) and the total gas energy (internal energy, e, kinetic energy, 
and magnetic energy) respectively, t is the time, p = pmnrtu is 
the mass density, fi = 1.3 is the mean atomic mass (assuming 
cosmic abundances), ra« is the mass of the hydrogen atom, «h 
is the hydrogen number density, u is the gas velocity, T is the 
temperature, and B is the magnetic field. We use the ideal gas 
law, P — (y - l)pe, where y — 5/3 is the adiabat ic index. The 



simul ations are performed using the flash code dFrvxell et al 



2000), an adaptive mesh refinement multiphysics code for as- 
trophysical plasmas. 

As initial cond i tions, we adopted the model profiles of 
Truelove & McKed |l999), assuming a spherical remnant with 



radius ro snr = 4 pc and with total energy Eg = 1.5 x 10 erg, 
originating from a progenitor star with mass of 15 M slm , 
and propagating through an unperturbed magnetohydrostatic 
medium. The initial total energy is partitioned so that 1/4 of 
the SN energy is contained in thermal energy, and the other 
3/4 in kinetic energy. The explosion is at the center (x, y, z) = 
(0, 0, 0) of the computational domain which extends between 
-30 and 30 pc in all directions. At the coarsest resolution, 
the adaptive mes h algorithm used in the flash code (paramesh; 
MacNeice et al . 2000) uniformly covers the 3-D computational 
domain with a mesh of 8 3 blocks, each with 8 3 cells. We al- 
low for 3 levels of refinement, with resolution increasing twice 
at each refine ment level. The refinement criterion adopted 
dLohnerll 1987b follows the changes in density and temperature. 
This grid configuration yields an effective resolution of w 0.1 
pc at the finest level, corresponding to an equivalent uniform 
mesh of 512 3 grid points. We assume zero-gradient conditions 
at all boundaries. 

We follow the expansion of the remnant for 22 kyrs, consid- 
ering two sets of simulations: 1) through a gradient of ambient 
density with a uniform ambient magnetic field; or 2) through 
a homogeneous isothermal medium with a gradient of ambi- 
ent magnetic field strength. Table Q] summarizes the physical 
parameters characterizing the simulations. 

In the first set of simulations, the ambient magnetic field 
is assumed uniform with strength B = 1 yuG and oriented 
parallel to the x axis. The ambient medium is modeled with 
an exponential density stratification along the x or the z direc- 
tion (i.e. parallel or perpendicular to the B field) of the form: 
n(£) = «o + «; exp(— f/ft) (where £ is, respectively, x or z) with 
no = 0.05 cm" 3 and n, = 0.2 ctrT 3 , and where h (set either to 
25 pc or to 10 pc) is th e density scale length. Th is configuration 
has been used by e.g. iHnatyk & Petrukl (119991) to describe the 
SNR expansion in an environment with a molecular cloud. Our 
choice leads to a density variation of a factor ~ 6 or ~ 60, re- 
spectively, along the x or the z direction over the spatial domain 
considered (60 pc in total). The temperature of the unperturbed 
ISM is T — 10 4 K at f = and is determined by pressure bal- 
ance elsewhere. The adopted values of T = 10 4 K, n — 0.25 
cm" 3 and B — 1 yuG at f = 0, outside the remnant, lead to 
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Table 1. Relevant initial parameters of the simulations: «o and 
rij are particle number densities of the stratified unperturbed 
ISM (see text), h is the density scale length, and (x,y,z) are 
the coordinates of the magnetic dipole moment. The ambient 
medium is either uniform or with an exponential density strat- 
ification along the x or the z direction (x-strat. and z-strat., 
respectively); the ambient magnetic field is uniform or dipolar 
with the dipole oriented along the x axis and located at (x, y, z). 



Model DZ1 



Model DX1 





ISM 


«o 


m 


h 


B 


(x, y, z) 






cirr 3 


cirr 3 


pc 




pc 


GZ1 


z-strat. 


0.05 


0.2 


25 


uniform 




GZ2 


z-strat. 


0.05 


0.2 


10 


uniform 




GX1 


x-strat. 


0.05 


0.2 


25 


uniform 




GX2 


x-strat. 


0.05 


0.2 


10 


uniform 




DZ1 


uniform 


0.25 






z-strat. 


(0,0,-100) 


DZ2 


uniform 


0.25 






z-strat. 


(0,0,-50) 


DX1 


uniform 


0.25 






x-strat. 


(-100, 0,0) 


DX2 


uniform 


0.25 






x-strat. 


(-50,0,0) 



B ~ 17 (where B = P/(B 2 /8n) is the ratio of thermal to mag- 
netic pressure) a ty pical order of magnitude of B in the diffuse 
regions of the ISM (Ma c Low & Klesserj|20o4 . 

In the second set of simulations, the unperturbed ambient 
medium is uniform with temperature T = 10 4 K and parti- 
cle number density n = 0.25 cirr 3 . The ambient magnetic 
field, B , is assumed to be dipolar. This idealized situation is 
adopted here mainly to ensure magnetostaticity of the non- 
uniform field. The dipole is oriented parallel to the x axis and 
located on the z axis (x = y = 0) either at z - -100 pc or at 
z = -50 pc; alternatively the dipole is located on the x axis 
(y — z = 0) either at x = -100 pc or at x = -50 pc (as shown 
in Fig. [TJ. In both configurations, the field strength varies by a 
factor ~ 6 (z or x = -100 pc) or ~ 60 (z or x — -50 pc) over 
60 pc: in the first case in the direction perpendicular to the av- 
erage ambient field (B), whereas in the second case parallel to 
(B). In all the cases, the initial magnetic field strength is set to 
B = 1 yuG at the center of the SN explosion (x — y — z = 0). 

Note that the trans ition time from adiabatic to rad iative 
phase for a SNR is (e.g. lBlondin et alJll998HPetrukll2005l) 



2.84 x 10 4 E 



4/17 n7 9 ' 11 yr . 

is:m J 



(5) 



where E$\ = Eo/(10 51 erg) and «j sm is the particle number den- 
sity of the ISM. In our set of simulations, runs GZ2 and GX2 
present the lowest values of the transition time, namely f tl * 25 
kyr. Since we follow the expansion of the remnant for 22 kyrs, 
our modeled SNRs are in the adiabatic phase. 

2.2. Nonthermal electrons and synchrotron emission 

We synthesize the radio emission from the remnant, assum- 
ing that it is only due to synchrotron radiation from relativistic 
electrons distributed with a power law spectrum N(E) = KE ^, 
where E is the electron energy, N(E) is the number of elec- 
trons per unit volume with arbitrary directions of motion and 
with energies in the interval [E, E + dE], K is the normaliza- 
tion of the electron distribution, and f is the power law index. 




10 10 

1 [pc] 



Figure 1. 2-D sections in the (x,z) plane of the initial mass 
density distribution and initial configuration of the unperturbed 
dipolar ambient magnetic field in two cases: dipole moment lo- 
cated on the z axis (DZ1, left panel), or on the x axis (DX1, 
right panel). The initial remnant is at the center of the domain. 
Black lines are magnetic field lines. 



Following lGinzburg & Syrovatskiil ( 1965 ). the radio emissivity 
can be expressed as: 



i(v) = dKBl +l v- a , 



(6) 



where C\ is a constant, B ± is the component of the magnetic 
field perpendicular to the line-of-sight (LoS), v is the frequency 
of the radiation, a = — 1)/2 is the synchrotron spectral index 
(assumed to be uniform everywhere and taken as 0.5 as ob- 
served in many BSNRs). To compute the total radio intensity 
(Stokes parameter I) at a given frequency vq, we integrate the 
emissivity /(vq) along the LoS: 



/(vo) 



z'(vo) dl , 



(7) 



where dl is the increment along the LoS. 

The normalization of the electron distribution K s in Eq. 
[6] (index "s" refers to the immediately post-shock values) de- 
pends on the injection efficiency (the fraction of electrons 
that move into the cosmic-ray pool). Unfortunately, it is un- 
known how the injection efficiency evolves in time. On theo- 
retical grounds, K s is expected to vary with the shock velocity 
V s h(f) and, in case of inhomogeneous ISM, with the immedi- 
ately post-shock value of mass d ensity, p s ' let us assume that 
approximately^ oc p^V^(t)~ b . iReynoldsl ( 1998 ) considered 
three empirica l alternatives for b as a fre e parameter, namely, 
b = 0. - 1 . -2. IPetruk & Bandieral d2006l) showed that one can 
expect b > and its value could be b m 4. Stronger shocks 
are more successful in accelerating particles. To be accelerated 
effectively, a particle should obtain in each Fermi cycle larger 
increase in momentum, which is proportional to the shock ve- 
locity. Negative b reflects an expectation that injection effi- 
ciency may behave in a way similar to acceleration efficiency: 
stronger shocks might inject particles more effectively. In con- 
trast, positive b represents a different point of view: efficien- 
cies of injection and acceleration may have opposite depen- 
dencies on the shock velocity. Stronger shock produces higher 
turbulence which is expected to prevent more thermal particles 
to recross the shock from downstream to upstream and to be, 
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therefore, injected. Since the picture of injection is quite un- 
clear from both theoretical and observational points of view, 
we do not pay attention to the physical motivations of the value 
of b. Instead, our goal is to see how different trends in evolu- 
tion of injection efficiency may affect the visible morphology 
of SNRs. Such understanding could be useful for future obser- 
vational tests on the value of b. 



We found, in agreement with iRevnoldsl (|1998), that the 
value of b does not affect the main features of the sur- 
face brightness distribution if SNR evolves in uniform ISM. 
Therefore we use the value b — to produce the SNR images in 
models with uniform ISM (models DZ1, DZ2, DX1, and DX2). 
In cases where non-uniformity of ISM causes variation of the 
shock velocity in SNR (models GZ1, GZ2, GX1, and GX2), 
we calculate im ages for b = -2, 0, 2. We follow the model of 
IRevnoldsl d 1 998b in description of the post-shock evolution of 
relativistic electrons. Adopting this approach and considering 
that ( = 2 (being a - 0.5, see above), one obtains that (see 
AppendixlAl 



K(a, t) 



K s {R,t) \P(R,t) 



P{a,t)\ lPa(ay 



p (R) 



p(a, t) \ 
p(R,t)j 



56/6+4/3 



(8) 



where a is the lagrangian coordinate, R is the shock radius, p is 
the gas density, P is the gas pressure, and the index "o" refers 
to the pre-shock values. It is important to note that this formula 
accounts for variation of injection efficiency caused by the non- 
uniformity of ISM. 

The electron injection efficiency may also vary with the 
obliquity angle between the external magnetic field and the 
shock normal, <pn n . The numerical simulations suggest that in- 
jection efficiency is larger for parallel shocks, i.e. where the 
magnetic field is parallel to the shoc k speed (obliquity an- 
gle close to zero; Elliso n et al. 1995). However, it has been 
shown dFulbright & Revnoldsl 1 19901) that models with injec- 
tion strongly favoring parallel shocks produce SNR maps that 
do not resemble any known objects (it is also claimed that 
injection is more efficient wher e the magnet ic field is per- 
pendicular to the shock speed; IJokipiil 119871) . On the other 
hand, comparison of known SNRs morphologies with model 
SNR images calculated for different strengths of the injec- 
tion efficiency dependence on obliquity suggests that the in- 
jection efficiency in real SNRs could not depend on obliquity 
(Petruk, in preparation). In such an unclear situation, we con- 
sider the three cases: quasi-parallel, qu asi-perpendicular, and 



isotro pic injection models. Following Fulbright & Reynolds 



ropi 

dl990l) . we model quasi-parallel injection by multiplying the 



nor malization of the ele ctron distribution K by cos 2 <p Bn 2 (see 



also Leckband etal . 1989), where <f> Bn 2 is the angle between the 
shock normal and the post-shock magnetic fielcQ. By analogy 
with the quasi-parallel case, we model quasi-perpendicular in- 
jection by multiplying K by sin 2 <fen2- 



1 For a shock compression ratio of 4 (the shock Mach number is » 
10 in all directions during the whole evolution in each of our simula- 
tions), the obliquity angle between the external magnetic field and the 



An important point is the degree of ordering of magnetic 
field downstream of the shock. Ra dio polarization ob servation 



of a number of SNRs (e.g Tvcho lDickel etal.ll 19911, SN1006 



shoc k normal, fan, is related to <ft Bn 2 by sin 
(e.g. lFulbright&Revnoldslll990h . 



(cot 2 Bn /i6+ir 



Revnolds & Gilmorelll993l) show the low degree of polariza- 
tion, 10-15% (in case of ordered magnetic field th e value ex- 
pected is about 70%: iFulbright & Reynoldlll990h . indicating 
highly disordered magnetic field. Thus we calculate the syn- 
chrotron images of SNR for two opposite cases. First, since our 
MHD code gives us the three components of magnetic field, 
we are able to calculate images with ordered magnetic field. 
Second, we introduce the procedure of the magnetic field dis- 
ordering (with randomly oriented magnetic field vector with 
the same magnitude in each point) and then synthesize the ra- 
dio maps. In models which have a disordered magnetic field, 
we use the post-shock magnetic field be fore disordering to cal- 
culate the angle (pBni', as discussed by IFulbright & Reynolds 
dl990h . this corresponds to assume that the disordering process 
takes place over a longer time-scale than the electron injection, 
occurring in the close proximity of the shock. Since we found 
that the asymmetries induced by gradients either of ambient 
plasma density or of ambient magnetic field strength are not 
significantly affected by the degree of ordering of the magnetic 
field downstream of the shock, in the following we will focus 
on the models with disordered magnetic field. 

The goal of this paper is to look whether non-uniform ISM 
or non-uniform magnetic field can produce asymmetries on 
BSNRs morphology. In order to clearly see the role of these 
two factors in determining the morphology of BSNRs, we use 
some simplifying assumptions about electron kinetic and be- 
havior of magnetic field in vicinity of the shock front. Our 
calculations are performed in the test-particle limit, i.e. they 
ignores the energy in cosmic rays. In particular, we do not con- 
sider possible amplifi cation of magnetic field by the cosmic-ra y 
streaming instability dLucek & BehfeOOdlBeil & LucekkOOll) . 
We expect that the main features of the modeled SNR morphol- 
ogy will not change if this process is independent of obliquity 
angle. If future investigations show undoubtedly that magnetic 
field amplification varies strongly with obliquity, the role of 
this effect in producing BSNRs have to be additionally studied. 

3. Results 

In all the models examined, we found the typical evolution 
of adiabatic SNRs expanding through an organized ambient 
magnetic field (see lBalsara et alJ l2001 and references therein): 
the fast expansion of the shock front with temperatures of few 
millions degrees, and the development of Richtmyer-Meshkov 
(R-M) instability, as the forward and revers e shocks progres s 
through the ISM and ejecta, respectively (see Ka ne et al . 1999). 
As examples, Fig. [2] shows 2-D sections in the (x, z) plane of 
the distributions of mass density and of magnetic field strength 
for the models GZ2, DZ2, and DX2 at t = 18 kyrs. The in- 
ner shell is dominated by the R-M instability that causes the 
plasma mixing and the magnetic field amplification. In the in- 
ner shell, the magnetic field shows a turbulent structure with 
preferentially radial components around the R-M fingers (see 
Fig. [3]). Note that, some authors have invoked the R-M insta- 
bilities to explain the dominant radial magnetic field observed 



S. Orlando et aL: On the origin of asymmetries in bilateral supernova remnants 



5 



Log Mass Density [gm cm"^ 



Log B-FieH [G] 




Log Mass Density [gm cm ] 



-10 10 
XltKl 



10 10 



Figure 2. 2-D sections in the (x, z) plane of the mass density 
distribution (left panels), in log scale, and of the distribution 
of the magnetic-field strength (right panels), in log scale, in 
the simulations GZ2 (upper panels), DZ2 (middle panels), and 
DX2 (lower panels) at t — 18 kyrs. The box in the upper left 
panel marks the region shown in Fig. [3] 



in the inner shell of SNRs (e.g. ljun & Normanll 996): however, 
in our simulations, the radial tendency is observed well inside 
the remnant and not immediately behind the shock as inferred 
from observations. 

We found that, throughout the expansion, the shape of the 
remnant is not appreciably distorted by the ambient magnetic 
field because, for the values of explosion energy and ambi- 
ent field strength (typical of SNRs) used in our simulations, 
the kinetic energy of the shock is many orders of magnitude 
larger than the energy dens ity in the ambient B field (see also 



Mineshige & Shibata Il990h ." The shape of the remnant does 



not differ visually from a sphere also in the cases with den- 
sity stratification of the ambient medium^ as it is shown by 
Hnatyk&Petrukl(ll999l) . 



The radio emission of the evolved remnants is character- 
ized by an incomplete shell morphology when the viewing an- 
gle is not aligned with the direction of the average ambient 



2 In these cases, the remnant appears shift ed toward the low den- 
sity re gion; see upper panels in Fig.[2](see also lDohm-Palmer & Jonesl 
fl996h . 




24.0 



24.5 



25.0 



25.5 



Figure 3. Close-up view of the region marked with a box in 
Fig. |2] The dark fingers mark the R-M instability. The mag- 
netic field is described by the superimposed arrows the length 
of which is proportional to the magnitude of the field vector. 



magnetic field (cf. lFulbright & Revnoldslll990l) : in general, the 
radio emission shows an axis of symmetry with low levels of 
emission along it, and t wo bright limbs (arcs) on either side 
(see also lGaenslerlll998l) . This morphology is very similar to 
that observed in BSNRs. 

3.1. Obliquity angle dependence 

For each of the models listed in Table Q] we synthesized the 
synchrotron radio emission, considering each of the three cases 
of variation of electron injection efficiency with shock obliq- 
uity: quasi-parallel, quasi-perpendicular, and isotropic particle 
injection. As an example, Fig. [4] shows the synchrotron radio 
emission synthesized from the uniform ISM model DZ1 with 
randomized internal magnetic field at t — 18 kyrs in each of 
the three cases. We recall that for these uniform density cases, 
we have adopted an injection efficiency independent from the 
shock speed (b = 0, Sect. 12.2b . All images are maps of total in- 
tensity normalized to the maximum intensity of each map and 
have a resolution of 400 beams per remnant diameter (Dsnr)- 
The images are derived when the LoS is parallel to the average 
direction of the unperturbed ambient magnetic field (B) (LoS 
aligned with the x axis), or perpendicular both to (B) and to 
the gradient of field strength (LoS along y), or parallel to the 
gradient of field strength (LoS along z). 

The different particle injection models produce images that 
can differ considerably in appearance. In particular, the quasi- 
parallel case leads to morphologies of the remnant not repro- 
duced by the other two cases: a center-brightened SNR when 
the LoS is aligned with x (top left panel in Fig. |4]i, a BSNR 
with two bright arcs slanted and converging on the side where 
B field strength is higher when the LoS is along y (top center 
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Figure 4. Synchrotron radio emission (normalized to the maximum of each panel), at t — 18 kyrs, synthesized from model DZ1 
assuming b — (see text) and randomized internal magnetic field, when the LoS is aligned with the x (left), y (center), or z (right) 
axis. The figure shows the quasi-parallel (top), isotropic (middle), and quasi-perpendicular (bottom) particle injection cases. The 
color scale is linear and is given by the bar on the right. The directions of the average unperturbed ambient magnetic field, (B), 
and of the magnetic field strength gradient, V|B |, are shown in the upper left and lower right corners of each panel, respectively. 



panel), and a remnant with two symmetric bright spots located 
between the center and the border of the remnant when the LoS 
is along z (top right panel). Neither the center-brightened rem- 
nant nor the double peak structure, showing no structure de- 
scribable as a shell, seems to be observed in We found 

analogous morphologies in all the models listed in TableQ] con- 
sidering the quasi-parallel c ase. As extensively discussed by 
Fulbright^Reynoldsl {l990) for models with uniform ambient 
magnetic field and b - -2, we also conclude that the quasi- 
parallel case leads to radio images unlike any observed SNR 
(see also lKesteven & Caswelllll987l) . 



The isotropic case leads to remnant's morphologies simi- 
lar to those produced in the quasi-perpendicular case although 
the latter case shows deeper minima in the radio emission than 
the first one. When the LoS is aligned with x (middle left and 
bottom left panels in Hig. SJ» or with y (middle center and bot- 
tom center panels), the remnants have one bright arc on the side 
where the B strength is higher. When the LoS is aligned with z 
(middle right and bottom right panels), the remnants have two 
opposed arcs that appear perfectly symmetric. We found that 
the isotropic and quasi-perpendicular cases lead to morpholo- 
gies of the remnants similar to those observed. 



3 Excluding filled center and composite SNRs, but these are due to 
energy input from a central pulsar. 
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Figure 5. Presentation as in Fig.|4]for model GZ1 with randomized internal magnetic field, assuming quasi-perpendicular particle 
injection and b — -2 (top panels), b — (middle) and b — 2 (bottom). The directions of the average unperturbed ambient magnetic 
field, (B), and of the ambient plasma density gradient, Vp, are shown in the upper left and lower right corners of each panel, 
respectively. 



3.2. Non-uniform ISM: dependence from parameter b 

For models describing the SNR expansion through a non- 
uniform ISM (models GZ1, GZ2, GX1, GX2), we derived the 
synthetic radio maps considering three alternatives for the de- 
pendence of the injection efficiency on the shock speed, namely 
b = -2,0,2 (see Sect. I2.21 i. As an example, Fig. [3] shows the 
synthetic maps derived from model GZ1 with randomized in- 
ternal magnetic field, assuming quasi-perpendicular particle in- 
jection, and considering b = -2 (top panels), b = (middle) 
and b - 2 (bottom). 

When the LoS is not aligned with the density gradient, the 
radio images show asymmetric morphologies of the remnants. 
In this case, the main effect of varying b is to change the de- 
gree of asymmetry observed in the radio maps. In the example 



shown in Fig.|3J the density gradient is aligned with the z axis 
and asymmetric morphologies are produced when the LoS is 
aligned with x (left panels) or with y (center panels). In all the 
cases, the remnant is brighter where the mass density is higher. 
On the other hand, the degree of asymmetry increases with in- 
creasing value of b. 

The reason of such behavior consists in the balance be- 
tween the roles of the shock velocity and of density in chang- 
ing the injection efficiency. Consider, as an example, the top 
left panel in Fig. [5] the increase of the shock velocity on the 
north (due to fall of the ambient density) leads to an increase of 
the brightness there (due to rise of the injection efficiency) that 
partially balances the increase of the brightness on the south 
due to higher density of ISM. On the other hand, for the model 
shown in the bottom left panel in Fig. [5] the fraction of accel- 
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Figure 6. Synchrotron radio emission (normalized to the maximum of each panel), at t — 18 kyrs, synthesized from models 
assuming a gradient of ambient plasma density (panels A and D; with b — 2) or of ambient magnetic field strength (panels B 
and E; with b — 0) when the LoS is aligned with the y axis. All the models assume quasi-perpendicular particle injection. The 
directions of the average unperturbed ambient magnetic field, (B), and of the plasma density or magnetic field strength gradient, 
are shown in the up per left and lower right co rner s of each panel, respectively. The right panels show two examples of radio maps 
(data adapted from Whi teoak & Green 1996 and iGaenslerl l 1998: the arrows point in the north direction) collected for the SNRs 
G338. 1+0.4 (panel C) and G296.5+10.0 (panel F). The color scale is linear and is given by the bar on the right. 



erated electrons increases on the south due to both the rise of 
density and the decrease of the shock velocity. 

When the LoS is aligned with the density gradient, the radio 
images are symmetric. In the example shown in Fig. |3J this 
corresponds to the maps derived when the LoS is along z (right 
panels); the remnants are characterized by two opposed arcs 
with identical surface brightness. 

3.3. Morphology 

Fig.|6]shows the radio emission maps, at a time of 18 kyrs, syn- 
thesized from models with a gradient of ambient plasma den- 
sity (panels A and D; assuming b — 2) and of ambient B field 
strength (panels B and E; assuming b = 0). All the models as- 
sume quasi-perpendicular particle injection (the isotropic case 
produces radio maps with similar morphologies and the quasi- 
parallel case is discussed later) and randomized internal mag- 
netic field. The viewing angle is perpendicular both to the av- 
erage direction of the unperturbed ambient magnetic field, (B), 
(direct along the x axis) and to the gradients of density or field 



strength (direct either along z, panels A and B, on, panels D 
and E). The right panels show, as example s, the radio maps of 
the S NRs G338.1+04 (panel C, data f romlWhiteoak & Green 
Il996l) and G296.5 + 10.0 (panel F, from lGaenslerlll998l) . 

In the quasi-perpendicular case discussed here, the max- 
imum synchrotron emissivity is reached where the magnetic 
field is strongly compressed. Th is configuration has bee n re- 
ferred as "equatorial belt" (e.g. iRothenflug et alj|2004 : (B) 
runs between the two opposed arcs (along the x axis). We found 
that, when the density or the magnetic field strength gradient is 
perpendicular to the field itself, the morphology of the radio 
map strongly depends on the viewing angle. In these cases, the 
two opposed arcs appear perfectly symmetric when the LoS 
is aligned with the gradient (see, for instance, the right pan- 
els in Fig. 0, otherwise the two arcs can have very different 
radio brightness, leading to strongly asymmetric BSNRs (see 
panels A and B in Fig. |6|. In the former case (LoS aligned 
with the gradient), the remnant is characterized by two axes of 
symmetry: one between the two symmetric arcs and the other 
perpendicular to the two. In models with strong magnetic field 
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strength gradients (DZ2; B varies by a factor ~ 60 over 60 pc), 
we found that the radio images are center-brightened when the 
LoS is aligned with the gradient (figure not reported). The fact 
that center-brightened remnants are not observed suggests that 
the external B varies moderately in the neighborhood of the 
remnants. 

In case of asymmetry, the gradient is always perpendicular 
to the arcs, and the brightest arc is located where either mag- 
netic field strength or plasma density is higher (see panels A 
and B in Fig. ©, since the synchrotron emission depends on 
the plasma density, on the pressure, and on the field strength 
(see Eqs. |6]and|8j>; in this case, there is only one axis of sym- 
metry oriented along the density or B gradient. When the LoS 
is parallel to (B) (along x in our models), the radio maps show 
a shell structure with a maximum intensity located where mag- 
netic field strength or plasma density is higher (see left pan- 
els in Fig. [4] for isotropic and quasi-perpendicular cases and 
left panels in Fig. |5j. Our simulations show that, when the 
density or the magnetic field strength gradient is perpendicu- 
lar to the field itself, remnants with a monopolar morphology 
can be observed at LoS not aligned with the gradient (see also 
Reynolds & Fulbrightll99oT) . Examples of observed monopolar 
remnants are G338. 1 +0.4 (see panel C in Fig.© or G327.4+ 1 .0 
orG341.9-0.3. 

When the density or B field strength gradient is parallel to 
(B) (panels D and E in Fig. [6j and the LoS lies in the plane per- 
pendicular to (B), the morphology of the radio map does not 
depend on the viewing angle and the two opposed arcs have 
the same radio brightness. In these cases, however, there is only 
one axis of symmetry and the two arcs appear slanted and con- 
verging on the side where field strength or plasma density is 
higher; again, the symmetry axis is aligned with the density 
or B strength gradient. Examples of this kind of objects are 
G296.5+10.0 (see panel F in Fig.© or G332.4-004or SN1006 
(which is, however, much younger than the simulated SNRs). 
When the external magnetic field is parallel to the LoS, because 
the system is symmetric about the magnetic field, the remnant 
is axially symmetric and the radio maps show a complete radio 
shell at constant intensity. 

In the quasi-parallel case, (B) runs across the arcs. This 
configuration has been referred as "polar caps" and it has been 
invoked for the SN1006 remnant dRofhenflug et al.ll2004 . The 
quasi-parallel case, apart from the center-brightened morphol- 
ogy discussed in Sect. 13. II can also produce remnant morpholo- 
gies similar to those shown in Fig. [6] As examples, Fig. [7] 
shows the radio emission maps obtained in the cases dis- 
cussed in Fig. ©but assuming quasi-parallel instead of quasi- 
perpendicular particle injection. Again, the viewing angle is 
perpendicular both to (B) (direct along the x axis) and to the 
gradients of density or field strength (direct either along z, pan- 
els A and B, or x, panels C and D). In the quasi-parallel case, 
remnants with a bright radio limb are produced if the gradi- 
ent of ambient density or of ambient B field strength is par- 
allel to (B) (instead of perpendicular to (B) as in the quasi- 
perpendicular case), whereas slanting similar radio arcs are ob- 
tained if the gradient is perpendicular to (B) (instead of parallel 
as in the quasi-perpendicular case). 



4. Discussion 

Our simulations show that asymmetric BSNRs are explained 
if the ambient medium is characterized by gradients either of 
density or of ambient magnetic field strength: the two opposed 
arcs have different surface brightness if the gradient runs across 
the arcs (see panels A and B in Fig. [6] and panels C and D in 
Fig. |7), whereas the two arcs appear slanted and converging 
on one side if the gradient runs between them (see panels D 
and E in Fig. ©and panels A and B in Fig. [7J. In all the cases 
(including the three alternatives for the particle injection), the 
symmetry axis of the remnant is always aligned with the gradi- 
ent. 

From the radio maps, we derived the azimuthal intensity 
profiles: we first find the point on the map where the intensity 
is maximum; then the contour of points at the same distance 
from the center of the remnant as the point of maximum in- 
tensity defines the azimuthal radio intensity profile. Following 
Fulbright & Reynolds! £1990), we quantify the degree of "bipo- 
larity" of the remnants by using the so-called azimuthal inten- 
sity ratio A, i.e. the ratio of maximum to minimum intensity 
derived from the azimuthal intensity profiles. In addition, we 
quantify the degree of asymmetry of the BSNRs by using a 
measure we call the azimuthal intensity ratio R m!lx > 1, i.e. 
the ratio of the maxima of intensity of the two limbs as de- 
rived from the azimuthal intensity profiles, and the azimuthal 
distance 8^, i.e. the distance in deg of the two maxima. In the 
case of symme tric BSNRs, R max = 1 and O n = 180°. As al- 
ready noted by Fulbright & Reyn olds ( 1990h . the parameter A 
depends on the spatial resolution of the radio maps and on the 
aspect angle (i.e. the angle between the LoS and the unper- 
turbed magnetic field); moreover we note that, in real observa- 
tions, the measure of A gives a lower limit to its real value if the 
background is not accurately taken into account. On the other 
hand, the parameters R mdx and 0q have a much less critical de- 
pendency on these factors and, therefore, they may provide a 
more robust diagnostic in the comparison between models and 
observations. 

Fig. [8] shows the values of A, R mdx , and 6>d derived for all 
the cases examined in this paper, considering the LoS aligned 
with the y axis, and radio maps with a resolution of 25 beams 
per remnant diametei0 (Dsnr)- Note that, our choice of the 
LoS aligned with y (aspect angle <p = 90°) implies that the 
values of A in Fig. [8] are upper limits, being A maximum at 
4> — 90" and minimum at <f> = 0° (see iFulbright & Reynolds 
1990h . The three models of particle injection (isotropic, quasi- 
perpendicular and quasi-parallel) lead to different values of A. 
In the isotropic and quasi-perpendicular cases, most of the val- 
ues of A range between 5 and 20 (for model DX2, A is even 
larger than 100); in the quasi-parallel case, the values of A are 
larger than 500. 

We found that, in general, a gradient of the ambient mag- 
netic field strength leads to remnant morphologies similar to 
those induced by a gradient of plasma density (compare, for 
instance, panel A with B and panel D with E in Fig. ©. On the 
other hand, if b < in GX and GZ models, ambient B field 



4 After the radio maps are calculated, they are convolved with a 
gaussian function with <r corresponding to the required resolution. 
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Figure 7. Presentation as in Fig. [6] as- 
suming quasi-parallel instead of quasi- 
perpendicular particle injection. 



gradients are more effective in determining the morphology of 
asymmetric BSNRs. This is seen in a more quantitative form 
in Fig. [8] DX and DZ models give R max values higher and 6 D 
values lower than GX and GZ models with b < 0: a modest 
gradient of the magnetic field (models DX1 and DZ1) gives 
a value of R max higher or 6>d lower than the two models with 
strong density gradients (models GX2 and GZ2) and b < 0. 

Fig. [8] also shows that, in models with a density gradient, 
the degree of asymmetry of the remnant increases with increas- 
ing value of b; the GX and GZ models with b > give values 
of R m -dx and 0£, comparable with (or, in the case of R mdx , even 
larger of) those derived from DX and DZ models. In the case of 
quasi-parallel particle injection for remnants with converging 
similar arcs, it is necessary a strong gradient of density perpen- 
dicular to B and b > (compare models GZ1 and GZ2 in the 
lower panel in Fig. [8]) to give values of 0£, comparable to those 
obtained with a moderate gradient of ambient B field strength 
perpendicular to B (see model DZ1 in Fig. [8]). 

In order to compare our model predictions with observa- 
tions of real BSNRs, we have selected 1 1 SNR shells which 
show one or two clear lobes of emission in archive total in- 
tensity radio images, separated by a region of minima. We 
have discarded all those cases in which several point-like or 
extended sources appear superimposed to the bright limbs, or 
other cases in which the location of maximum or minimum 
emission around the shell is difficult to derive. Unlike other lists 



1987tlFuibright & Revnoldsll990tlGaenslerlll 998). here we fo- 
cus on a reliable measure of the parameters A, R max and 6^; we 
avoid, therefore , patchy and irre gular limbs, as in the case of 
G320.4-01.2 of iGaenslerl dl998l) . Moreover, we are obviously 
not dis carding remnants which have constraints on A, R m!ix or 
#d (e.g. Fulbright & Reynolds! 1990l considered only cases with 
Rmax < 2), and we are considering remnants observed with a 
resolution greater than 10 beams per remnant diameter. Since 
in our models we follow the remnant evolution during the adi- 
abatic phase, we also need to discard objects that are clearly 
in the radiative phase. Unfortunately, for most of the objects 
selected, there is no indication of their evolutionary stage in lit- 
erature. Assuming that the remnant expands in a medium with 
particle number density n; sm < 0.3 crrT 3 , the shock radius de- 
rived from the Sedov solution at time t a (i.e. at the transition 
time from the adiabatic to the radiative phase; see Eq. [5]) is 



r tt = 19 E 



5/17 -7/17 



5] 



n 



< 35 pc , 



(9) 



of BSNRs published in the literature (e.g. lKesteven & Casw ell 



where we have assumed that £51 = 1.5. Therefore, we only 
considered remnants with radius r snr < 35 pc (i.e. with size < 
70 pc) that are, most likely, in the adiabatic phase. Our list does 
not pretend to be complete or representative of the class, and 
it is compiled to derive the observed values of the parameters 
A, R mdx and 0£, with the lowest uncertainties. For this reason, 
we have considered remnants for which a total intensity radio 
image in digital format is available. Actually, in most of the 
cases, we have use d the 843 MHz data of the MOST supernova 
remnant catalogue (Whit eoak & Greenlll996l) . 
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Figure 8. Azimuthal intensity ratio A (i.e. 
the ratio of maximum to minimum inten- 
sity around the shell of emission - see text; 
upper panel), azimuthal intensity ratio /f max 
(i.e. the ratio of the maxima of intensity 
of the two limbs around the shell; middle 
panel), and azimuthal distance 9 D (i.e. the 
distance in deg of the two maxima of in- 
tensity around the shell; lower panel) for 
all the cases examined, considering the LoS 
aligned with the y axis and a spatial res- 
olution of 25 beams per remnant diame- 
ter, £>snr- Black crosses: isotropic; red tri- 
angles: quasi-perpendicular; blue diamonds: 
quasi-parallel. 



Our list is reported in Table [2] We have separated evolved 
and young SNRs. While the young SNRs listed in Table|2]have 
very reliable measurement of A, R mdx and 0£, and a good record 
of literature, making them very good candidate to test the diag- 
nostic power of our model, we stress that the models we are 
considering in this paper are focused on evolved SNRs; we 
leave the discussion about young SNRs to a separate work. For 
each object in Table [2] we show the apparent size, the distance 
(from dedicated studies where possible, otherwise fro m the re- 
vised S — D relation of Case & B hattacharva 1998; see their 
paper for caveats on usage of the 2 - D relation to derive SNR 
distance), the real size, the resolution of the observation, and 
the parameters A, R msix , and 8o we have introduced here. 

Table |2] shows that most of the 11 remnants have A < 10, 
i.e. values consistent with those derived in Fig. [8]for the three 
alternatives for the particle injection (recall that the values 
shown in the figure have to be considered as upper limits). 
Four remnants show high values ofA(10<A< 100) that 
are difficult to explain in terms of the isotropic or the quasi- 
perpendicular injection models with b < unless the remnant 
expands through a non-uniform ambient magnetic field (see 
models DX2, and DZ2 in Fig. [8]). In the light of these consid- 
erations, we cannot exclude a priori any of the three alternative 
models for the particle injection. 



Four of the 1 1 objects in Table [2] show values of R mdx > 2, 
pointing out that, in these objects, the bipolar morphology is 
asymmetric with the two radio limbs differing significantly in 
intensity. An example of this kind of remnants is G338. 1+0.4 
(see panel C in Fig. [6j. In the light of our results, its morphol- 
ogy can be explained if a gradient of ambient density or of 
ambient magnetic field strength is either perpendicular to the 
average ambient magnetic field, (B), in the isotropic and quasi- 
perpendicular cases or parallel to (B) in the quasi-parallel case. 
It is worth noting that reveling such a gradient from the obser- 
vations may be a powerful diagnostic to discriminate among 
the alternative particle injection models, producing real ad- 
vances in the understanding of the nonthermal physics of strong 
shock waves. 

An extreme example of a monopolar remnant with a bright 
radio limb is G327.4+1.0 whose value of R mllx is larger than 
10. Fig. [8] shows that high values of R max can be easily ex- 
plained as due to non-uniform ambient magnetic field strength 
or to non-uniform ambient density if b > 0. We suggest that the 
morphology of G327.4+1.0 may give some hints on the value 
of b (and, therefore, on the dependence of the injection effi- 
ciency on the shock velocity) if the observations show that the 
asymmetry is due to a non-uniform ambient medium through 
which the remnant expands. 
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Table 2. List of barrel-shaped SNR shells for which a measurement of A, R m . dx , and 6 D is presented for comparison with our 
models. 



Remnant^ 


Flux 


size 


d 


size 


Res.^ 


A 




e D 


Ref./notes 




Jy 


arcmin 


kpc 


pc 


beams/DsNR 






deg 




Evolved Remnants 


G296.5+10.0 


48 


90x65 


2.1 


55 x40 


108 


> 11 


1.2 


85 


1 


G299.6-0.5 


1.1 


13 x 13 


18.1 


68 


18 


6 


2 


160 


2 


G304.6+0.1 


18 


8x8 


7.9 


18 


11 


20 


1.5 


120 


3 


G327.4+1.0 


2.1 


14x13 


13.9 


56 


19 


> 10 


> 10 


ND 


2,4 


G332.0+0.2 


8.9 


12 x 12 


< 20 


< 70 


17 


5 


1 


145 


2,7 


G338.1+0.4 


3.8 


16 x 14 


9.9 


46x40 


21 


3 


2 


> 120 


2 


G341. 9-0.3 


2.7 


7x7 


14.0 


28 


10 


8 


3 


170 


2 


G346.6-0.2 


8.7 


11 x 10 


8.2 


26x23 


15 


2 


1.1 


110 


2.7 


G351.7+0.8 


11 


18 x 14 


6.7 


35 x27 


22 


2 


1.6 


130 


2 










Young 


Remnants 










G327.6+14.6 


19 


30x30 


2.2 


19 x 19 


42 


22 


1 


127 


5 


G332.4-0.4 


34 


11 x 10 


3.1 


10x9 


15 


7 


1.6 


98 


6 



Reference s and notes. - (1) A.k . a. PK S 1209-51/52. Age: 3-20 kvrs. iRoger et al.l dl988h- Distance fromlGiacani et alj d200ol) . (2) Distance 
derived bv lCase & Bhattacharval i 19981) usi ng a revised X - D relation. (3 ) Distance from Caswell et al.l dl975» . (4) This shell has only one limb 
("monopolar" according to the defin ition of lFulbright & Revnoldsll 19901). A and R mdx are lowe r limits and no do is derived. (5) A.k.a. SN1006. 
Distance from I Winkler et al.l J2003h . (6) A.k.a. RCW103. Distance from lRevnoso et alJ ( l2004T) . (7) Two maxima have been found in one lobe. 
#d is the average of the two. 



" All the data are from the MOST supernova remnant catalogue jWhiteoak & GreerJll996h . except where noted. 
* Spatial resolution of the observation in beams per remnant diameter. 



In Table|2] six of the 1 1 remnants (including the two young 
remnants SN1006 and RCW103) have values of 6 D < 140°, 
pointing out that, in these objects, the two bright radio limbs 
appear slanted and converging on one side. An example of this 
class of objects is G296.5+10.0 (a.k.a PKS 1209-51/52) shown 
in panel F in Fig. [6] In this case, the value of 6t> ~ 85° de- 
rived from the observations may be easily explained as due 
to a gradient of magnetic field strength either parallel to (B) 
in the isotropic and quasi-perpendicular cases or perpendicu- 
lar to (B) in the quasi-parallel case. Models with a gradient 
of ambient density cannot explain the low values of 8 D found 
for G296. 5+10.0 unless the gradients are strong (the density 
should change by a factor 60 over 60 pc) and the dependence 
of the injection efficiency on the shock velocity give^J/? > 2. 

5. Conclusions 

Our findings have significant implications on the diagnostics 
and lead to several useful conclusions: 

1. The three different particle injection models (namely, 
quasi-parallel, quasi-perpendicular and isotropic dependence 
of injection efficiency from shock obliquity) can produce con- 
siderably different images (see Fig. |4j». The isotropic and quasi- 
perpendicular cases lead to radio images similar to those ob- 



5 Large positive values of b do not necessarily mean an increas- 
ing fraction of shock energy going into relativistic particles as the 
shock slows down because decelerating shock accelerates particles to 
smaller £ max , namely the maximum energy at which the electrons are 
accelerated. 



served. The parallel-case may produce radio images unlike any 
observed SNR (center-brightened or with a double-peak struc- 
ture not de scribable as a shell). This is in agreement with the 
findings of Fulbri ght & Reynolds 



2. In models with gradients of the ambient density, the 
dependence of the injection efficiency on the shock velocity 
(through the parameter b defined in Sect. 12. 2t affects the degree 
of asymmetry of the radio images: the asymmetry increases 
with increasing value of b. 

3. Small variations of the ambient magnetic field lead to 
significant asymmetries in the morphology of BSNRs (see 
Figs. |6] and 17). Therefore, we conclude that the close similar- 
ity of the radio brightness of the opposed limbs of a BSNR 
(i.e. R mdx ~ 1 and 8u ~ 180") is evidence of uniform ambient 
B field where the remnant expands. 

4. Variations of the ambient density lead to asymmetries 
of the remnant with extent comparable to that caused by non- 
uniform ambient magnetic field if b = 2. 

5. Strongly asymmetric BSNRs (i.e. R max » 1 or #d ^ 
180°) imply either moderate variations of B or strong (moder- 
ate) variations of the ISM density if b < 2 (b > 2) as in the 
case, for instance, of interaction with a giant molecular cloud. 

6. BSNRs with different intensities of the emission of the 
radio arcs (i.e. R mdx > 1) can be produced by models with a 
gradient of density or of magnetic field strength perpendicular 
to the arc (upper panels in Fig. [6] and lower panels in Fig. |7j, 
and the brightest arc is in the region of higher plasma density 
or higher magnetic field strength. 
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7. Remnants with two slanting similar arcs (i.e. 6 D < 180°) 
can be produced by models with a gradient of density or of 
magnetic field strength running centered between the two arcs 
(lower panels in Fig. [6] and upper panels in Fig. [7]), and the 
region of convergence is where either the plasma density or the 
magnetic field strength is higher. 

8. In all the cases examined, the symmetry axis of the rem- 
nant is always aligned with the gradient of density or of mag- 
netic field. 

We found that the degree of ordering of the magnetic field 
downstream of the shock does not affect significantly the asym- 
metries induced by gradients either of ambient plasma density 
or of ambient magnetic field strength; thus our conclusions, de- 
rived in the case of disordered magnetic field, do not change in 
the case of ordered magnetic field. 

We defined useful model parameters to quantify the degree 
of asymmetry of the remnants. These parameters may provide a 
powerful diagnostic in the comparison between models and ob- 
servations, as we have shown in a few cases drawn from a ran- 
domly selected sample of BSNRs presented in Table [2] For in- 
stance, if the density of the external medium is known by other 
means (e.g. thermal X-rays, HI and Co maps, etc.), BSNRs can 
be very useful to investigate the variation of the efficiency of 
electron injection with the angle between the shock normal and 
the ambient magnetic field or to investigate the dependence of 
the injection efficiency from the shock velocity. Alternatively, 
BSNRs can be used as probes to trace the local configuration 
of the galactic magnetic field if the dependence of the injection 
efficiency from the obliquity is known. 

It is worth emphasizing that our model follows the evolu- 
tion of the remnant during the adiabatic phase and, therefore, 
their applicability is limited to this evolutionary stage. In the 
radiative phase, the high degree of compression suggested by 
radiative shocks leads to increase of the radio brightness due 
to compression of ambient magnetic field and electrons. Since 
our model neglects the radiative cooling it is limited to rela- 
tively small compression ratios and, therefore, it is not able to 
simulate this mechanism of limb brightening. 

It will be interesting to expand the present study, consider- 
ing the detailed comparison of model results with observations. 
This may lead to a major advance in the study of interactions 
between the magnetized ISM and the whole galactic SNR pop- 
ulation (not only BSNRs), since the mechanisms at work in the 
BSNRs are also valid for SNRs of more complex morphology. 
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Appendix A: Derivation of Eq. © 

We follow Reynolds (1998) in the description of the evolution of elec- 
tron distribution. His approach is extended here to the possibility to 



deal with non- uniform ISM (cf. Petru kll2006h . Fluid element a = R(ti) 
was shocked at time t,, where R is the radius of the shock, and a is the 
Lagrangian coordinate. At that time the electron distribution on the 
shock was 



(A.l) 



where E\ is the electron energy at time r ; , K s is the normalization of 
the electron distribution immediately after the shock (in the following, 
index "s" refers to the immediately post-shock values), and f is the 
power law index. Since we are interested in radio emission, we have 
to account for only energy losses of electrons due to the adiabatic 
expansion jRevnoldsll 19981) : 



dE E dp 
dt 3p dt ' 

where p is the mass density, so the energy varies as 



E = Ei 



P(a,t) 



1/3 



(A.2) 



(A.3) 



The conservation law for the number of particles per unit volume per 
unit energy interval 



N(E,a,t) = N(Ei,a,ti) 



a 2 da dE; 
err 2 dr dE ' 



(A.4) 



where cr is the shock compression ratio and r is the Eulerian coor- 
dinate, together with the continuity equation p (a)a 2 da = p(a, t)r 2 dr 
(index "o" refers to the pre-shock values) and the derivative 



dE, I p(a, t) 



1/3 



dE \Ps(a, ti) j 
implies that downstream 



N(E,a,t) = K s (a,h)E 



/ \ \ 1 



(f+2)/3 



(A.5) 



(A.6) 



If A" s oc Pf,V sh (t) h , where V s h(f) is the shock velocity and p s is the 
immediately post-shock value of density, then 



Therefore, the distribution of relativistic electrons follows 
K(a, t) N(E, a, t) ! Po (a) \ I V sh (t) \ b j p(a, t) * ( ' : ' 



K s (R,t) N(E, R, t) \p (R)l\V sh ( tl )l \p s (a,td 



(A.7) 



(A.8) 



Now we can substitute Eq, IA.8l w ith the ratio of the shoc k velocities 
which comes from the expression dHnatvk & Petrukll 19991) 



P(a,t) _ /PcW r 2/3 ( V&(ti) \ 2 l p(a,t) ^ 1 



P s (R,t) \Po(R)) \V sh (t)j \p s (R,t) 
Thus, finally 



F( ., I D t A\-W 2 / ( N\-(&+f-0/3 / , A \5/>/6+Gr+2)/3 

K(a,t) I r(a,t) \ I p (a) \ / p{a,t) 1 



K s (R,t) \P S (RJ)I \Po(R) 



P S (RJ) 



(A.9) 



(A. 10) 



This formula may easily be used to calculate the profile of K(a) for 
known P(a) and p(a) in the case of the radial flow of fluid. In the case 
when mixing is allowed, the position R should correspond to the same 
part of the shock which was at a at time t\. 
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